Effect of bulk charged impurities on the bulk and surface transport in 
three-dimensional topological insulators 

Dedicated to the memory of Professor Anatoly Larkin 

Brian Skinner, Tianran Chen, and B. I. Shklovskii 
Fine Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 55455, USA 

(Dated: April 1, 2013) 

In the three-dimensional topological insulator (TI), the physics of doped semiconductors exists 
literally side-by-side with the physics of ultra-relativistic massless Dirac fermions. This unusual 
pairing creates a novel playground for studying the interplay between disorder and electronic trans- 
port. In this mini-review we focus on the disorder caused by the three-dimensionally distributed 
charged impurities that are ubiquitous in TIs, and we outline the effects it has on both the bulk and 
surface transport in TIs. We present self-consistent theories for Coulomb screening both in the bulk 
and at the surface, discuss the magnitude of the disorder potential in each case, and present results 
for the conductivity. In the bulk, where the band gap leads to thermally activated transport, we 
show how disorder leads to a smaller-than-expected activation energy that gives way to VRH at low 
temperatures. We confirm this enhanced conductivity with numerical simulations that also allow 
us to explore different degrees of impurity compensation. For the surface, where the TI has gapless 
Dirac modes, we present a theory of disorder and screening of deep impurities, and we calculate 
the corresponding zero-temperature conductivity. We also comment on the growth of the disorder 
potential as one moves from the surface of the TI into the bulk. Finally, we discuss how the presence 
of a gap at the Dirac point, introduced by some source of time-reversal symmetry breaking, affects 
the disorder potential at the surface and the mid-gap density of states. 



I. INTRODUCTION 

The three-dimensional (3D) topological insulator 
(TT)pmi has generated a great deal of excitement in the 
physics community because of its gapless surface states, 
wh ich host a spectrum of quantum transport phenom- 
engPU. Unfortunately, while a number of crystals have 
been identified to be 3D TIs, most of them are not ac- 
tually insulators, but instead have a relatively large bulk 
conductivity that shunts the surface conductivity for TI 
crystals of substantial thickness (> 10 /mi). How to 
achieve a bulk-insulating state is a problem that is widely 
discussed in the current literature^E^. 

Typically, as-grown TI crystals are heavily doped re- 
type semiconductors, so that the Fermi level resides in 
the bulk conduction band. In order to arrive at a bulk 
insulating state, such TIs are compensated by acceptors. 
With increasing compensation K — Na/Nd, where Nr> 
and Na are the concentrations of monovalent donors and 
acceptors, respectively, the Fermi level shifts from the 
conduction band to inside the gap and then into the va- 
lence band. When compensation of donors is complete, 
K = 1, the Fermi level is in the middle of the gap and 
the most insulating state of the TI is reached. The hope 
is that for a TI with bulk band gap E g ~ 0.3 eV (as, for 
example, in Bi2Se3) the bulk resistivity should obey the 
activation law 

p = po exp(A/fc B T) (1) 

with activation energy A = E g /2 ~ 0.15 eV, so that at 
room temperatures and below the TI is well insulating. 

The typical experimental situation near K = 1, how- 
ever, is frustrating^. In the range of temperatures be- 
tween 100 K and 300 K the resistivity is activated, but 



with an activation energy that is roughly three times 
smaller than expected, A ~ 50meV. At T ~ 100 K the 
activated transport is replaced by variable range hop- 
ping (VRH) and the resistivity grows even more slowly 
with decreasing T. Finally, at even smaller temperature, 
T < 50 K, the resistivity saturates^ at a value < 10 ilem. 

In a recent papei^ we showed that the unexpectedly 
large bulk conductivity of TIs at K = 1 can be explained 
as a consequence of the enormously-fluctuating random 
Coulomb potential created by randomly-positioned donor 
and acceptor impurities. In later papers we extended this 
analysis to the case of near complete compensation^ 
K < 1 and 1 — K <C 1, and we examined the effect of 
random Coulomb impurities on the surface disorder and 
transport properties^. In this mini-review our goal is 
to outline in a general way the effects of random, 3D- 
distributed Coulomb impurities in TIs on both the bulk 
and surface properties. We describe the screening mech- 
anisms for the random Coulomb potential both within 
the bulk of the TI and at the surface, and we present 
predictions for the magnitude of the disorder potential 
and the conductivity. 

Our theoretical treatment is also motivated by the re- 
cent experiments of Ref. 2 1 where the random potential 
at the surface of typical TIs (Bi2Sc3 and Bi2Te3) was 
studied directly by spectroscopic mapping with a scan- 
ning tunneling microscope. It was shown that near the 
Dirac energy random fluctuations of the potential have 
a Gaussian-like distribution with a width ~ 20 - 40 meV 
that can be attributed to deep impurity charges. We 
show below that such fluctuations are consistent with dis- 
order produced by three-dimensionally distributed bulk 
Coulomb impurities that are screened by the gapless TI 
surface. 
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Crucial to our theoretical description throughout this 
paper is the assumption of a random spatial distribu- 
tion of impurities. This assumption is readily justified 
for TI samples made by cooling from a melt, where the 
distribution of impurities in space is a snapshot of the 
distribution that impurities have at higher temperature, 
when their diffusion practically freezes^. In 3D TIs, as 
in conventional narrow band gap semiconductors, at this 
temperature there is a concentration of intrinsic carri- 
ers larger than the concentration of impurities. Intrinsic 
carriers thus screen the Coulomb interaction between im- 
purities, so that impurities remain randomly distributed 
in space. When the temperature is lowered to the point 
where intrinsic carri ers re combine, the impurities are left 
in random position d 23 l 24 [ If the diffusion of impurities 
freezes at T ~ 1000 K it is reasonable to assume that im- 
purities are randomly positioned for semiconductors with 
bulk band gap E g < 0.3 eV. Throughout this paper we 
deal with such narrow band gap TIs, such as Bi2Se3, for 
which our description of randomly-positioned impurities 
is accurate. We also assume everywhere that donor and 
acceptor energy levels are shallow, meaning that their 
binding energy is much smaller than E g . 



The remainder of this paper can be divided into two 
parts. In the first part, comprising Sees. |H|- |IV[ we fo- 
cus on bulk properties, essentially treating the TI as a 
strongly- or completely-compensated semiconductor and 
ignoring the surface states. In Sec. |n]we give a concep- 
tual explanation of the bulk disorder potential and the 
origin of the anomalously small bulk resistivity. Sec. |III| 
formulates a numerical model of the TI bulk and uses it 
to calculate the corresponding electron density of states 
(DOS). In Sec. IV we present our algorithm for comput- 



ing the thermally activated conductivity, analyze our re- 
sults, and arrive at an expression for the unusually small 
bulk activation energy. We also evaluate the localization 
length of states close to the Fermi energy and estimate 
the characteristic temperature associated with variable- 
range hopping. 



The second part of this paper, comprising Sees. [V] - 
|VIII[ deals with the effects of Coulomb impurities on the 
properties of the TI surface. In Sec. [V] we describe a self- 
consistent theory of the screened disorder potential at 
the TI surface and compare it with experiment. Sec. |VII| 
uses this theory to calculate the conductivity of surface 
electrons. Sec. |VI| briefly discusses how the amplitude 
of the disorder potential transitions from its large bulk 
value to its smaller value at the surface. Finally, Sec. 
|VIII| discusses an extension of our analysis to the case 
where the TI surface has a gap introduced by some source 
of time-reversal symmetry breaking. Where applicable, 
the major results of each section are summarized at the 
beginning of the section. 



II. ORIGIN OF THE ENHANCED BULK 
CONDUCTIVITY 

As mentioned in the Introduction, randomly- 
positioned impurities create a disordered Coulomb 
landscape in the bulk of the TI, which has the effect 
of reducing the activation energy A relative to what 
one would naively expect by thinking about flat valence 
and conduction bands. In this section we explain this 
idea more fully, focusing first on the case of complete 
compensation, where the bulk transport can be de- 
scribed using the th eory of a completely compensated 
semiconductor (CCS^MH. 

This theory is based on the idea that at K — 1, when 
almost all donors and acceptors are charged, random spa- 
tial fluctuations of the local concentration of impurities 
result in large fluctuations of charge. Their potential is 
poorly screened, because of the vanishing average concen- 
tration n = Nd~Na of electrons, and therefore has huge 
fluctuations. These fluctuations bend the conduction and 
valence band edges and in some places bring them to the 
Fermi level, creating electron and hole puddles that in 
turn non-linearly screen the random potential. As a re- 
sult, the amplitude of potential fluctuations is limited by 
Eg/2, so that the ground state, illustrated sch ema tically 
in Fig. [TJ resembles a network of p-n junction d 24 * 25 [ The 
characteristic size of these p-n junctions is^ 
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which can be thought of as the correlation length of 
the random potential. For the typical parameters E g ~ 
0.3 eV, Nd — 10 19 cm~ 3 , and dielectric constant K = 30, 
this length scale R « 150 nm > N D 1/3 = 4.6 nm. That 
is, we deal with a very long range potential. 

As a result of these long range fluctuations, the resis- 
tivity can be dramatically different from the naive ex- 
pectation based on thinking about flat valence and con- 
duction bands. First, at relatively high temperatures ac- 
tivated conductivity is due to electrons and holes acti- 
vated from the Fermi level to their corresponding classi- 
cal percolation levels (classical mobility edges), E e and 
Eh, in the conduction and the valence bands. According 
to numerical modeling^ at K = 1, the activation energy 
A ~ 0.15E g , meaning that E e and E^ are substantially 
closer to the Fermi level \i than to the unperturbed bot- 
tom of the conduction band, E c , or ceiling of the valence 
band, E v (Fig. [T^,). (E c and E v are the energies of 
the conduction and valence bands, respectively, as they 
would be in the absence of a random potential.) Thus, 
one can think of the universal small factor A/E g w 0.15 
as corresponding to a percolation threshold associated 
with percolation through the potential created by ran- 
dom Coulomb impurities in 3D. 

Second, at sufficiently low temperatures electrons and 
holes can hop (tunnel) directly between puddles, so that 
activated transport is replaced by VRH. In Ref. [T5] we 
showed that with decreasing temperature the activated 
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resistivity crosses over directly to the Efros-Shklovskii 
(ES) law^ 



p = Po exp(T ES /T) 



1/2 



(3) 



where Xes = Ce 2 /fcsK£, e is the electron charge, £ is the 
localization length of electron states with energy close to 
the Fermi level, and C w 4.4 is a numerical coefficient. 
Together our results for the activated and VRH resis- 
tivity established the universal upper limit of the bulk 
resistivity p(T) for a 3D TI compensated by shallow im- 
purities. 
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FIG. 1. Energy diagram of a) a completely compensated semi- 
conductor (K — 1) and b) a strongly compensated semicon- 
ductor (1 — K <C 1) with band gap E g . The upper and the 
lower straight lines indicate the unperturbed positions of the 
bottom of the conduction band, E c , and the ceiling of the 
valence band, E v ; the middle straight line corresponds to the 
Fermi level p. Meandering lines represent the band edges, 
which are modulated by the fluctuating potential of charged 
impurities. R is the characteristic size of potential fluctu- 
ations. Percolation levels (mobility edges) for electrons, E e 
and holes, Eh are shown by dashed lines. Puddles occupied by 
carriers are shaded. Shallow impurities levels are not shown 
because they practically merge with band edges. 



K = 1 exactly, and it is important to know how the 
results for a CCS, where K = 1, are extended to the 
case of a strongly compensated semiconductor (SCS), for 
which < 1 — K -C 1. For example, one can ask at 
which value of 1 — K does the activation energy A be- 
come twice smaller than at K = 1. For definiteness we 
consider an n-type SCS, where the concentration of elec- 
trons n — Nd — Na <C Ne,. We model numerically the 
ground state of such a SCS and its resistivity using al- 
gorithms similar to those of Ref. US We find that, in 
agreement with analytic theory^, when 1 — K grows the 
screening of the random potential improves and the cor- 
relation length R of the random potential decreases. The 
amplitude of the random potential decreases as well; hole 
puddles shrink and eventually vanish; and the chemical 
potential p moves up, so that E c — p decreases. One can 
say that with increasing (1 — K) screening happens by 
bending of the conduction band only, while all acceptors 
remain occupied by electrons and negatively charged. All 
these changes are illustrated by the transition from a) to 
b) in Fig. [l] 

As a result of these changes with growing 1 — K, the 
activation energy A decreases. We find that the relation 
A = 0.3(E C — p) obtained in Ref. [TS]for K = 1 remains 
valid for 1 — K <C 1 as well (see Fig. [6] below) . [In 
p-type semiconductors, where K — Nd/Na, & similar 
relationship holds: A = 0.3(p - E v ).] By K = 0.97 
the activation energy A is already several times smaller 
than at K = 1. This result shows that achieving the 
maximum bulk resistivity, with A = 0.l5E g , is not easy. 
It also helps to explain the origin of the large scatter in 
the magnitude of A among different TI samples^. 

Our prediction that A = 0.3(E C — p) can in principle 
be directly compared with experiments in TIs. Indeed, 
for each K the position of the Fermi level, (E c — p), 
can be found via measurements of the concentration of 
electrons in the surface states using Shubnikov-dc-Haas 
oscillations. 

At lower temperatures the activated bulk conduction 
crosses over to ES VRH. In Sec. IV we study this 
crossover numerically and also show how Tes, which is 
correlated with A, decreases with I — K . 

It should be mentioned that these results for the bulk 
conductivity are also applicable to other narrow gap 
semiconductors, for example, to InSb. Historically, a 
large effort was made to make InSb insulating via strong 
compensation, with the goal of improving the perfor- 
mance of InSb-based photodetectors. Results were again 
frustrating: the dark resistivity was too small. Our re- 
sults are in reasonable agreement with transport experi- 
mental data for InStP^. 



In a Ref. 1191 we expand our focus to consider not just 
the maximum possible bulk resistivity that appears at 
K = 1, but to address the more practical question of 
the dependence of the bulk resistivity on the degree of 
compensation K at 1 — K <C 1. Indeed, with exist- 
ing methods of growth of TI samples one cannot get 



III. MODEL OF BULK IMPURITIES AND THE 
DENSITY OF STATES 

In order to study numerically the bulk properties of 
a heavily doped SCS, we introduce a model of the bulk 
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donors and acceptors. In this section we first describe our 
numerical model and then use it to calculate the position 
of the Fermi level relative to the band edges as a function 
of compensation, K, and to evaluate the density of states 
of impurity states. Our major results are shown below in 
Figs. |] and || 

Specifically, we model the bulk as a cube containing 
a large number of randomly-positioned donors and ac- 
ceptors. We numerate all donors and acceptors by the 
index i and use n, = or 1 to denote the number of 
electrons residing on a donor or acceptor. We also in- 
troduce the binary variable fi to discriminate between 
donors (for which /j = 1) and acceptors (fi = —1). The 
Hamiltonian of our system is then 

H = E -y/^i + E v ( r ij)m, (4) 

* (ij) 

where = (fi/2 — m + 1/2) is the net charge of site i, 
V(r) is the interaction energy between two like-charged 
impurities at a distance r, and all energies are defined 
relative to the Fermi level. The first term of Eq. Q con- 
tains the difference between the energies of donors and 
acceptors, which for the case of shallow impurities is very 
close to the semiconductor gap E g . The second term of 
H represents the total interaction energy of charged im- 
purities. Note that Eq. Q does not include the kinetic 
energy of electrons and holes in the conduction and va- 
lence bands and, therefore, aims only at a description of 
the low temperature physics of SCS (ksT <C E c — /i). 

The form of the interaction law V(r) requires some 
consideration. For two impurities at a distance r > ob, 
where as is the effective Bohr radius of impurity states, 
one can use for V(r) the normal Coulomb interaction 
V(r) = e 2 /nr. For example, one can consider a pair 
of empty and distant donors. In such a donor pair one 
donor shifts the energy of the electron level on the other 
by V(r) = — e 2 /nr. This classical form for V(r) is good 
for a lightly doped SCS. In a heavily doped SCS, on the 

— 1/3 

other hand, where as > N D , most impurities have 
at least one neighbor at distance r < a^, and quantum 
mechanical averaging over the electron wave function be- 
comes important. (This is why an uncompensated heav- 
ily doped semiconductor is a good metal.) For example, 
a pair of donors cannot create an electron energy state 
deeper than that of the helium-like ion, which has bind- 
ing energy 2e 2 /kcib is the binding energy of the shallow 
donor state. The interaction law V(r) should therefore be 
"softened" at short distances r < as to reflect quantum 
mechanical effects. We model this behavior by contin- 
uing to use the classical Hamiltonian of Eq. Q with a 
truncated Coulomb potential V(r) = e 2 / ' k{t 2 + a 2 ^) 1 / 2 . 

Below it is convenient to express energies in units of 
e 2 Np 3 /n. In these units, a typical TI with band gap 
0.3 eV has E g s» 30. We unfortunately could not model 
E g = 30 directly, since in this case the very large correla- 
tion length of the random potential, R, leads to large 
size effects. Instead, we present results for the more 



modest value E g = 15, for which the size effect requires 
extrapolation^ only for K = 1. Results for the smaller 
E g = 10 are largely identical. 

In our numerical simulations, we first randomly place 
donors and acceptors within the simulation volume; re- 
sults presented below correspond to M — 20000 donors 
and 200007^ acceptors. We then search for the arrange- 
ment of electrons (or equivalcntly, the set of electron oc- 
cupation numbers {n-i}) that minimizes H, and we use 
this set to calculate the DOS and the conductivity. We 
begin our search from the state where all MK acceptors 
are populated by electrons and negative (m = 1, (ft = 
— 1), and where an equal number of randomly chosen 
donors are empty and positive (rij = 0, (ft = 1), while 
the remaining M(l — K) donors are filled and neutral 
(rii = l,qi = 0). The charged donors and acceptors in 
this initial state create a random potential whose magni- 
tude exceeds E g , and as a result the system's energy is 
well above that of the ground state. In order to bring the 
system closer to its ground state, we attempt sequentially 
to transfer electrons from an occupied impurity (either a 
neutral donor or a negatively charged acceptor) to an 
unoccupied one (a positively charged donor or a neutral 
acceptor) . If the proposed move lowers the total system 
energy H, then it is accepted, otherwise it is rejected. To 
check whether H goes down with each proposed move, for 
a given set of electron occupation numbers {rii} it is con- 
venient to introduce the single-electron energy state, £j, 
at a given impurity i: 

, '[p. E ri 'V- ( 5 ) 

In the ground state, single electron energies must satisfy 
the ES criterion 

ej -Si- V{r l3 ) > (6) 

for all i, j with n,- = 1 and rij = 0. We use our numeri- 
cal simulation to loop through all pairs of impurity sites 
i, j and enforce this criterion; if a given pair does not 
satisfy Eq. ^ , then we move the electron from impurity 
i to j and recalculate all e.j. This process is continued 
until no single-electron transfers are possible that lower 
H . The final arrangement of electrons can be called a 
pseudo-ground state, since higher order stability crite- 
ria of the true ground state (corresponding to simultane- 
ously changing three or more electron numbers) are not 
checked. Such pseudo-ground states are known to de- 
scribe the propert ies of real ground states with a high 
degree of accurac} l 24 * 29 [ Results below are obtained at 
Eg = 15 and a B = N D 1/3 for K = 1,0.99,0.98,0.97,0.96 
and 0.95, and are averaged over 100 realizations of the 
impurity coordinates. 

For each pseudo-ground state we estimate the Fermi 
energy fi as the arithmetic average of the minimum empty 
and maximum occupied energies e. The results are shown 
in Fig. [2j which shows how the Fermi level n(K) shifts 
from the middle of the gap toward the conduction band 
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FIG. 2. (Color online) Distance between the Fermi level \i 
and the bottom of the conduction band E c as a function of 
1 — K , as calculated by numerical simulation. Energies are in 
units of e 2 Np 1 ^ 3 /k, and the simulated band gap is E g — 15. 
The size of dots characterizes the numerical uncertainty. 
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FIG. 3. (Color online) Dimensionless single-electron DOS 
g(e), in units of [(1 + K)N/ (e 2 N 1/3 / k)], as a function of elec- 
tron energy e calculated from the Fermi level. Results are 
plotted for K = 0.95 (blue) and K = 1 (red) using £ 9 = 15. 
Impurity states with e < are occupied and those with e > 
are empty. At K = 1 the total DOS of impurities has donor- 
acceptor symmetry, which is lost with growing 1 — K. 



bottom with growing 1 — K. At 1 — K > 0.01 this depen- 
dence is in reasonable agreement with the prediction of 
single band theory (which ignores the valence band and 
acceptors}^ that E c — fi = A(l — if) -1 / 3 , where A is 
a numerical coefficient. Note, however that for heavily 
doped SCS the coefficient Ah ~ 1.4 is twice smaller than 
the coefficient Ai ~ 2.8 obtained in Ref. [5Hfor a lightly 
doped SCS, for which Nna% <C 1. In the latter case the 
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short range Coulomb interaction at distances r -C N D 
leads to an additional contribution to n of the same order 
of magnitude. 

The resulting DOS of impurities is shown in Fig. [3] for 
K = 1 and K = 0.95. g*(e) is the DOS in the units of 
(1 + K)Nd I '(e 2 Np 3 1 ' k) and is normalized to unity. At 
K = 1, the nearly constant and symmetric DOS between 
£ = —Eg and e = E g reflects the practically uniform dis- 
tribution of the random potential from —E g /2 to E g /2 
and, correspondingly, of the band edges E c and E v be- 
tween to Eg and between to — E g) respectively (see 
Fig. [iji). Near the Fermi level (e = 0) one can see the 
ES Coulomb gar 1251 



On the other hand, at K < 1 the DOS of impurity 
states loses the donor-acceptor symmetry it has at K = 1. 
As described in Sec. [n] (see Fig. [IJ, with growing 1 — 
K hole puddles are eliminated so that acceptors become 
disengaged from screening. The acceptor DOS (leftmost 
peak of Fig. [3]) therefore splits from the donor one, which 
in turn develops two peaks separated by the Fermi level 
at e = 0. The large right peak belongs to empty donors, 
while the small and narrow left peak belongs to occupied 
donors (electron puddles). These two donor peaks are 
separated by the ES Coulomb gap. 



IV. NUMERICAL MODELING OF 
THERMALLY ACTIVATED CONDUCTIVITY 

In the previous section we described our procedure for 
finding the energy levels of donor and acceptor impuri- 
ties in the pseudo-ground state. We now discuss how 
these results can be used to calculate the bulk conduc- 
tivity of a SCS, and we present results for the conduc- 
tivity both in the high-temperature, activated regime 
and in the low-temperature, VRH regime. Our ma- 
jor results are twofold. First, we find that in the ac- 
tivated regime the activation energy decreases as the 
chemical potential approaches the conduction band ac- 
cording to A fa 0.3(E C — /x) [see Fig. [fjj. Second, wc 
study how the characteristic temperature Tes in the 
VRH regime depends on compensation, and we find that 

T ES ~ 4.4 v /A( e 2A^ /3 /K). 

Our process for numerically calculating the resistivity 
is as follows. Once the energies {e^} are known (as cal- 
culated using the procedure described in Sec. Ill), we 



evaluate the resistivity usi ng th e approach of the Miller- 
Abrahams resistor networlP^QJ. In this description each 
pair of impurities i,j is said to be connected by a link 



with resistance Rij = i?oexp[27'i 



the activation energy e,-j is defined^ as follows: 

\Sj - £i\ - V(rij), £j£i < 



Sij/ksT], where 



max £, 



EjEi > 0. 



(7) 



The resistivity of the system as a whole is found us- 
ing a percolation approach^!. Specifically, we find the 
minimum resistance R c such that if all links with resis- 
tance > R c are cut, then there still exists a percola- 
tion pathway connecting opposite faces of the simulation 
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FIG. 4. (Color online) The temperature dependence of the 
resistivity in the whole temperature range 200 > T* > 
0.03. The dimensionless resistance (lnp)* is plotted against 
(T*)~ 1//2 to illustrate that the resistivity follows the ES law 
at low temperatures. The dashed lines are the best linear fits. 
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FIG. 5. (Color online) The temperature dependence of the 
resistivity in the high temperature range 200 > T* > 1. The 
dimensionless resistance (In p)* is plotted against (T*) -1 to il- 
lustrate that the resistivity is activated at high temperatures. 
The dashed lines are the best linear fits. 



volume. This approach captures the exponential depen- 
dence of the resistivity on the temperature, and we ig- 
nore details of the prefactor. Below we plot the temper- 
ature in the dimensionless units T* — 2/csTft/e 2 A r ^/ 3 £ 
and the resistivity p using the dimensionless quantity 

1 /3 

(lnp)* = (£,Np /2) \nR c /R . These dimensionless units 
eliminate any explicit dependence on the localization 
length £. 

In Fig.Hlthe resulting resistivity is plotted as a function 
of (j 1 *)-™ over t ne huge range of temperature 200 > 
T* > 0.03 for four different values of the compensation 
K . The resulting linear dependence at 0.3 > T* > 0.03 
indicates that at low temperatures the resistivity is well 



The higher tem- 
otted separately as 



described by the ES law [see Eq. (|3) 
perature range 200 > T* > 1 is pi 
a function of 1/T* in Fig. [5j Here the linear slope 
suggests a well-defined activation energy that depends 
on the compensation K. At extremely high T* > 50, 
which generally corresponds to unrealistically large tem- 
peratures, the conduction is dominated by activation of 
carriers across the band gap, which is not captured by 
our model. 

Extracting the slope of the curves in Fig. [5] (dashed 
lines) gives an estimate of the activation energy A as a 
function of compensation K. Combining this result with 
the values for the chemical potential n(K) calculated in 
Sec. |III| yields the data shown in Fig. [6) where A is plotted 
as a function of (E c — /i) for all the studied values of 
compensation K = 1,0.99,0.98,0.97,0.96,0.95. One can 
see that the equation A ~ 0.3(£' c — /i) holds reasonably 
well for all K in this interval. 

So far we have emphasized results that do not explic- 
itly depend on the localization length £. In fact, £ de- 
termines the magnitude of Tes , an d therefore determines 
the value of temperature at which the conduction tran- 
sitions from activated to VRH behavior. We argue now 




FIG. 6. (Color online) The activation energy A as a func- 
tion of the distance between the Fermi level and the conduc- 
tion band, plotted for K = 1.0,0.99,0.98,0.97,0.96, and 0.95 
(from right to left). The dashed line is the best linear fit, 
A ~ 0.3(i5 c — fi). All energies are plotted in units of the band 
gap E g . 



that in a TI £ is quite large, leading to a prominent role 
for VRH. To see this, one can imagine an electron with 
energy close to the Fermi level tunneling from one elec- 
tron puddle to another, distant one. If such an elec- 
tron were to tunnel along the straight line connecting 
the two puddles it would tunnel through high barriers 
and its wave function would decay sharply, with a de- 
cay length ^ < as. However, this straight line does 
not constitute the path of least action for the tunnel- 
ing electron. Instead, a tunneling electron can use the 
same geometrical path as a classical percolating electron, 
which has energy A above the Fermi level, and thereby 
avoid large barriers. One can roughly estimate the tun- 
neling decay length by assuming that along such a "per- 
colating" tunneling path the potential energy barriers V 
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are uniformly distributed in the range < V < A and 
neglecting the additional contribution to the action as- 
sociated with curvature of this path. Integration over 
V then gives a localization length £ ~ H/ymA and 
fesT ES = 4.4(mA) 1 / 2 (e 2 /K/i). For a TI with a B = N D 1/3 

this implies k B T ES = AA^J A(e 2 N^ /3 /k). 

The dependence Tes oc \fK implies that when A in- 
creases ~ 2.5 times, as in Fig. [6] corresponding to the 
difference between K = 0.95 and K = 1, the ES temper- 
ature Tes increases by ~ 60%. For a TI with k — 30 and 
Np = 10 19 cm -3 , this corresponds to a variation in Tes 
from 500 to 800 K. The regime of ES VRH in TIs can be 
studied experimentally, but such a study requires suffi- 
ciently thick samples that the bulk conduction provides 
a larger contribution to the total conductance than the 
TI surfaces. 



V. SELF-CONSISTENT THEORY OF THE 
SURFACE DISORDER POTENTIAL 

In the first part of this paper we showed how the bulk 
conduction is strongly influenced by the presence of ran- 
dom Coulomb impurities, which produce large bending of 
the bulk conduction and valence bands. We now turn our 
attention to the problem of how these same impurities af- 
fect the surface transport provided by the Dirac-like sur- 
face states. For this problem we adopt the same model of 
monovalent Coulomb impurities that are randomly dis- 
tributed throughout the bulk of the TI, and we focus our 
attention on the case of complete (or nearly-complete) 
compensation Njj = Na = N, where the Fermi level 
lies within the bulk band gap. As we show below, for 
determining the properties of the surface one can safely 
ignore the weak nonlinear screening by electron and hole 
puddles formed in the bulk (illustrated in Fig. [I]). 

In this section we present a self-consistent theory for 
the magnitude of the disorder potential at the TI surface, 
following Ref. 1201 Our primary result is an expression 
for the amplitude of fluctuations of the electric potential 
energy, T, at the TI surface as a function of the chemical 
potential, fi, measured relative to the Dirac point. In 
particular, for /i = we show below that 



holes) per unit area in surface puddles is given by 



r 2 = 



Q 4/3 



(M = 0). 



(8) 



Here a = e 2 / k s Hv is the effective fine structure constant, 
where k s is the effective dielectric constant at the surface 
and v is the Dirac velocity. This expression describes 
screening of the disorder potential via the formation of 
electron and hole puddles at the TI surface. The charac- 
teristic size of these puddles is given by 



AT-l/3 



(P = 0), 



(9) 



2 2 / 3 a 4 /3' 

and the corresponding total number of electrons (or 



(m = o). 



(10) 



Eqs. Q - ( 10 ) are derived below, along with results cor- 
responding to large u. Below we also derive a simple 
relation for the autocorrelation function of the potential 
at the TI surface, which has an unusually slow decay and 
can be used to verify the bulk origin of disorder. These 
results were confirmed by numerical simulation in Ref. 

EDI 

Our primary tool for describing screening of the elec- 
tric potential is the Thomas-Fermi (TF) approximation, 
which applies in the limit where the potential <j>(r) varies 
slowly compared to the characteristic Fermi wavelength 
of electrons at the surface. Specifically, the TF approxi- 
mation gives 



H = E f [n(r)} - e<j>(r), 



(11) 



where Ef(n) = 
(e 2 / an s ) \J 47r|rt| sgn(n) is the local Fermi energy 
and n(r) is the 2D electron concentration at the point 
r on the surface. The TF approximation is justified 
whenever a <C 1, as we show below. In TIs such small 
a can be seen as the result of the large bulk dielectric 
constant k > 30. We note here that for describing 
the properties of the surface state, which exists at a 
dielectric discontinuity, one should use for the effective 
dielectric constant k s the arithmetic mean of the internal 
and external dielectric constants. If the TI is in vacuum, 
then k s = (« + l)/2 ~ re/2. 

When the chemical potential is large enough in magni- 
tude that /i 2 3> e 2 (4> 2 ), where (...) denotes averaging over 
the TI surface, the relation Ef(n) can be linearized to 
read Ef[n(r)} ~ u+Sn(r)/v(fi). Here Sn(r) = n(r)—no is 
the difference in the electron concentration relative to the 
state with zero electric potential, no = a 2 k 2 /! 2 /(Aire 4 ), 
and v(u) = a 2 K 2 |/i|/(27re 4 ) is the density of states at 
Ef = u. From this density of states one can define 
a screening radius r s = K s /2ire 2 v = e 2 /a 2 K, s jjL that 
characterizes the distance over which fluctuations in the 
Coulomb potential are screened by the surface. The 
TF approximation is valid when the Fermi wavelength 
■\f ~ n o ~ e 2 /an s u is much smaller than r s , which 
gives the condition a«l. 

One can understand qualitatively the magnitude of the 
potential fluctuations, T, using the following simple ar- 
gument. For a given point on the TI surface, one can 
say that only impurities within a distance Bl < r s con- 
tribute to the potential; those impurities at a distance 
R' 3> r s are effectively screened out (one can say that 
they are screened by their image charges in the "metal- 
lic" TI surface). Impurities with R' < r s , on the other 
hand, are essentially unscreened. There are ~ iVr 3 such 
impurities, and their net charge is of order Q ~ e-JNr^, 
with a random sign. The absolute value of the potential 
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at the surface is then ~ Q/K s r s , so that T ~ eQ/n s r s ~ 
{e 2 N 1 / 3 / Ks )(Nr 3 s ) 1 ^ ~ y/e 2 N/ Ka v ~ y/e^/a 2 ^. 

In order to more accurately derive the value of T, one 
can start by considering the potential created by a single 
impurity charge +e. When such an impurity charge is 
placed a distance z from the TI surface (say, above the 
origin), it creates a potential <j>i(r; z) that within the TF 
approximation is given by^ 



exp[— qz\ 



k s J 1 + (qr s )- 



Y J (qr)dq, (12) 



where Jo ( x ) is the zeroth order Bessel function of the first 
kind. At large z/r s , Eq. (12) can be expanded to give 



<j)i(r;z) 



k s (r 2 + z 2 ) 3 / 2 



A simple physical derivation of Eq. ( 13 ) is based on the 
notiorP^l that for a distant impurity, such that z r s , a 
surface with screening radius r s effectively plays the role 
of a metallic surface positioned below the real surface at 
a distance z = — r s /2. Equation (13) can then be viewed 



formation of electron and hole puddles at the surface. 
Nonetheless, one can define an average density of states 
(v) at the surface, which determines, self-consistcntly, 
the typical screening radius r s and the magnitude of the 
potential fluctuations at the TI surface. This value (v) 
can be equated with the thermodynamic density of states 
of the system, d/j,/d(n), where (n) is the overall electron 
concentration of the surface. 

Consider, for example, the case fx = 0, where by sym- 
metry the average value of the potential {4>) — 0. At 
any given point r on the surface, the potential </>(r) is 
the sum of contributions from many individual impurity 
charges, provided that the characteristic screening radius 
r s = K s /2Tre 2 (v) > A -1 / 3 . This implies that, by the 
central limit theorem, the value of the potential across 
the surface is Gaussian-distributed with some variance 
(13) (0 2 ) = T 2 /e 2 that remains to be calculated. Within the 
TF approximation the local density of states at the point 
r is u[— e<j>{r)] = ea 2 K 2 |</>(r)|/(27re 4 ), so that one can cal- 
culate the average density of states as 



as the sum of the potentials created by the original charge 
at a distance z above the plane and its opposite image 
charge at a distance z + r s below the plane, expanded to 
lowest order in r s /z. 

The total potential at the origin is (f)(0) = 
<li4>i{ r ii z i)> where the index i labels all impurity 
charges, qt is the sign of impurity i, and and z% are the 
radial and azimuthal coordinates of its position. Under 
the assumption that all impurity positions are uncorre- 
cted and randomly-distributed throughout the bulk of 
the TI, the average of 4> 2 is given by 



l (r';z')] 2 2Nd 2 r'dz'. 



(14) 



Here, the quantity 2Nd 2 r'dz' describes the probability 
that the volume element d 2 r'dz' contains an impurity 
charge, and the integration is taken over the semi-infinite 
volume of the bulk of the TI. The width of the disorder 
potential at the TI surface, T, is defined by T 2 — e 2 ((f) 2 ). 
Inserting Eq. ( 12 ) into Eq. ( 14 ) and taking the integral 
then gives 



e 2 N 2ne 4 N 



K S V a 2 K 3 |/i|' 



Ml > 



e 2 N l/3 



K s a 



2/3 



(15) 



Eq. ([15]) is correct so long as the fluctuations in the 
Coulomb potential energy are small compared to the 
chemical potential, or T < \fj,\; this gives the condition 
written in parentheses. 

On the other hand, when is very small, the fluctu- 
ations in the Coulomb potential become large compared 
to the chemical potential, and one cannot talk about a 
spatially uniform local density of states v or screening 
radius r s . Instead, the Fermi energy has strong spatial 
variations, and the random potential is screened by the 




(16) 



This result for (v) can be inserted into the first equality 
of Eq. (15), T 2 = e 2 N/n s {v) 1 to give a self-consistent re- 
lation for the amplitude of potential fluctuations^. This 
procedure gives the result first announced at the be- 
ginning of this section, Eq. jSJ. Substituting Eqs. ^ 
and ( 16 1 into the expression for the screening radius, 
r s = K~/27re 2 (i/), gives Eq. 

One can also calculate the total concentration of elec- 
trons/holes in surface puddles, n p , implied by this result 
for r 2 . This is done by first inverting the TF relation, Eq. 

(a 2 K 2 /47re 2 ^ 2 - 



(11), at (j, = to give n(</>) = (a nj/Ane )</> sgn(^). In- 
tegrating this expression for n((j)) weighted by the Gaus- 
sian probability distribution for (f> gives 



exp [- 



Q 2 K 2 r 2 

87re 4 : 



( M = 0). 



/2r 2 

2 



Substituting the result of Eq. Q for T 2 then gives Eq. 
(10). One can also combine this result for the residual 
electron/hole concentration, n p , with the expression for 



the screening radius, r s , to arrive at an estimate for the 
number of electrons/holes per puddle: M p ~ im p r 2 ~ 
7r/16a 2 . Apparently at small a puddles typically contain 
many electrons/holes, M p ^> 1. 

Our primarily results, outlined in Eqs. ((8|) - (10), are 



valid within the TF approximation so long as the typical 
Fermi wavelength, A/ ~ e 2 /aK s T, is much smaller than 
the typical screening radius, r s ~ e 2 /a 2 K s r, which again 
gives the condition a«l. 
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As we mentioned above, at fx = the screening ra- 
dius r s describes the characteristic size of electron or hole 
puddles at the TI surface. More generally, r s plays the 
role of a length scale over which potential fluctuations 
at the surface are correlated. Such correlations can be 
discussed in a quantitative way by defining the potential 
auto-correlation function: 



C(r) = (0(RW + r)) t 



(17) 



where (...} r ' denotes averaging over the spatial coordi- 
nate r', and where by symmetry the correlation function 
depends on |r| — r only. In the remainder of this section 
we derive approximate analytical results for C(r), and 
show that spatial correlations in the potential have an 
unusually slow de cay. 

At r = 0, Eq. (17) reproduces the definition of (<ji> 2 ), 



scanning tunneling microscopy, one can discriminate be- 
tween disorder by bulk impurities and disorder by impu- 
rities located in a layer close to the surface. 

We now discuss the magnitude of T and r s implied by 
these expressions for typical TIs, which generally have an 
impurity concentration N ~ 10 19 cm" 3 . Typical values 
of the Dirac velocity and fine structure constant for TIs 
can be taken from Ref. [2TJ which reports Kv — 1.3 eV A 
and estimates a = 0.24. Using these parameters gives 
r ~ 30meV and r s ~ 20 nm at the Dirac point, /i = 0. 
At large \n\ > 30meV, both F 2 and r s decay as l/|/x|. 

As mentioned in the Introduction, the theory presented 
in this section provides a good description of the recent 
experimental results of Ref. [2"T1 where the random poten- 
tial at the surface of the 3D TIs Bi 2 Se 3 and Bi 2 Te 3 was 
studied using a scanning tunneling microscope^. Indeed, 
in these experiments it was found that the electric poten- 
tial at the surface was well-characterized by a Gaussian 
distribution with a standard deviation r ~ 10 - 20meV, 
and the characteristic length scale of potential fluctua- 
tions was estimated as r s ~ 20 - 30 nm. One can com- 
pare these measurements to our theoretical predictions 
by using the parameters listed above and inserting the 
measured chemical potential fj, ~ 100 meV into Eq. (15 1. 
This procedure gives T ~ 18 meV, and the corresponding 
screening radius r s ~ 5 nm, so that our theory is indeed 
in reasonably good agreement with experiment. Further, 
Ref. found that the disorder potential at the surface 
was not correlated with the position of surface impurities, 
indicating that the surface disorder potential originates 
primarily from impurities deep below the TI surface, as 
we have described. 

Throughout this section, we have worked within the 
assumption that bulk impurities are completely ionized, 
or in other words that there is no screening by conduction 
band electrons or valence band holes in the bulk. Such 
an assumption is valid when the chemical potential re- 
sides in the middle of a large bulk band gap. In this case 
donors or acceptors can only be neutralized by very large 
band bending discussed in Sec. |n] (see Fig.[l]). Such fluc- 
tuations take place over a long length scale R that scales 
as the square of the distance between the Fermi level and 
the nearest band edge [see Eq. ^] and is typically on the 
order of hundreds of nanometers for typical TliP^. On 
the other hand, near the surface of the TI the potential 
fluctuations are screened much more effectively and over 
a much shorter distance, r s , by the (ungapped) surface 
states. As shown above, r s is typically < 20 nm, and the 
amplitude of surface potential fluctuations r ~ 30 meV 
Eq. (fl9|) implies an unusually slow decay of potential < E g ~ 300 meV. One can therefore safely assume that 

near the surface there is no large band bending and one 
can indeed treat bulk impurities as completely ionized. 
The effect of bulk screening should appear only in the 
long-range behavior of the correlation function, r » JJ, 
where the 1/r decay of C(r) is truncated and, as one can 



so that C(0) = r ye 2 . At small enough distances that 
r « r s , one can expect that the value of C(r) is deter- 
mined primarily by unscreened impurities that are within 
a distance r s from the surface, as explained above during 
the derivation of T 2 . On the other hand, at r > r s cor- 
relations are produced primarily by impurities that are 
relatively far from the surface, as can be seen from the 
following scaling argument. Consider two surface points 
separated by a distance r»r s . One can imagine draw- 
ing a cube of size r that extends into the bulk of the TI 
and which contains the two surface points on opposite 
edges of one of its faces. Such a cube contains ~ Nr 3 im- 
purities , and has a net impurity charge with magnitude 
q ~ e\J Nr 3 and random sign. These impurity charges 
are located at a mean distance ~ r 3> r s above the sur- 



face and, therefore, by Eq. (13), contribute a net poten- 
tial ~ qr s / n s r 2 ~ (e/ ' n s )\J 'Nr 2 jr to both surface points. 



The square of this potential roughly gives the autocorre- 
lation of the potential, C{r) ~ e 2 Nr 2 / n 2 r . 

A more careful expression for C'(r) can be derived by 
writing 



C{r) = J 0i(r';z')0i(r' - r; z')2Nd 2 r' dz' , 



(18) 



similar to Eq. (14). Inserting the asymptotic expression 
of Eq. ( 13 1 for <\>\ and evaluating the integral gives 



, 2^e 2 A^ s 2 r 2 /e 2 

C r ~ ^ = r/r s >l). (19) 

K z r r/r s 

This result was also confirmed by numerical simulation 
in Ref. EH 



correlations at the surface, which, as explained above, 
arises from long-range fluctuations of the potential cre- 
ated by deep bulk impurities. This behavior can be 
contrasted with the much faster decay of C(r) that 
would result from a two-dimensional (2D) distribution of 
Coulomb impurities at a distance d from the surfaced 
C(r) ~ e 2 n i dr 2 /K 2 r' i , where n$ is the 2D impurity con- 
centration. Thus, by studying C(r) experimentally by 



show, is replaced with C(r) ~ e 2 NRr 2 /K 2 r 2 . 

Finally, we note that our theory ignores the possibility 
of screening by material outside the TI. For example, if 
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the TI is placed next to a metal electrode or an ionic 
liquid 3 -^, then this external material can screen the large 
potential fluctuations created by the bulk, thereby de- 
creasing r and r s . 



VI. FROM SURFACE TO BULK 

In Sec. [n] we showed that deep within the bulk of the 
TI the disorder potential has large fluctuations of order 
r ~ E„. On the other hand, in Sec.fVlwe showed that at 



. In 1 

the TI surface the disorder potentialhas a much smaller 
amplitude, T ~ (e 2 N 1 ^ 3 / na 2 ^ 3 ). In this section we elab- 
orate briefly on the crossover between these two results, 
or in other words we describe how the amplitude of po- 
tential fluctuations grow as one moves from the surface 
of the TI into the bulk. 

Generally speaking, as one moves a distance z > 
into the bulk of the TI, the amplitude of the disorder 
potential increases in magnitude. In order to see quan- 
titatively how r grows as a function of z, one can as- 
sume, for the moment, that the TI surface is equivalent 
to a perfect metallic plane. In this case, each impurity 
at position (r',z') has a corresponding image charge at 
(r', —z'), and the total potential at (0, z) is equal to the 
sum of the potentials created by the original impurity 
and its image. One can calculate T 2 (z) by averaging the 
square of this potential over all possible positions of the 



impurity charge [as in Eq. (14)]. This calculation gives 
r 2 (z) = 87riVe 4 z/K 2 . That is, T 2 (z) grows linearly with 
the distance z from the TI surface. This growth contin- 
ues until r becomes large enough that T 2 (z) = (E g /2) 2 , 
at which point electron and hole puddles begin to form in 
the bulk and one arrives at the bulk screening picture de- 
scribed in Ref.lTS"! This distance corresponds to z — i?/4; 
at smaller z the potential fluctuations are small enough 
that practically all donors and acceptors are charged. 

One can now recall that the TI surface is not perfectly 
metallic, and that its screening length r s is finite, so that 
r 2 (z) should be somewhat larger. In fact, at z ^> r s one 
can still use the formula above for T 2 (z) by introducing a 
small modification allowing for the fact that the metallic 
surface is effectively shifted to the position z = —r s /2 
(as discussed in Sec. [V]). Making this adjustment gives 



r 2 (z) = 87rA^e 4 (z + rJjTj/n 2 at z > r s , which does not 
significantly alter our conclusions. 



VII. SURFACE CONDUCTIVITY 

We now turn our attention to the problem of how the 
3D-distributcd Coulomb impurities within the TI bulk 
affect the surface conductivity. As discussed at the be- 
ginning of Sec. |V| we limit our consideration to the case 
where the Fermi level resides within the bulk band gap, 
where one can safely assume that all relevant bulk impu- 
rities are ionized. 



Our primary result is an expression for the electron 
conductivity a of the surface as a function of the average 
2D surface electron concentration n. In particular, when 
p is the typical puddle concentration 



n 3> n p , where n 



at /j, = [see Eq. (10)], we find that the conductivity is 
given by 



2\/i 



,3/2 



h a 2 ln(l/a) N 



(20) 



where e 2 /h is the conductance quantum. At much 
smaller electron concentrations, n <C n p , the conductiv- 
ity saturates at a value er m j n , which we estimate as 



1 



h ira ln(l/a) 



(21) 



To derive these results, we first note that in the limit of 
large chemical potential fi, where the electron density is 
only weakly modulated by the disorder potential, one can 
show using the Boltzmann kinetic equation that for elec- 
trons with a massless Dirac spectrum the conductivity is 
given bySsHSU 



e 2 /it 



(22) 



Here r is the momentum relaxation time. In the limit of 
zero temperature, the scattering rate 1/t can be found by 
integrating the squared scattering potential produced by 
a given impurity over all impurities and over all scattering 
angles. More simply, one can arrive at an expression 
for 1/t by taking the result for the scattering rate of a 
2D layer of impurities with concentration m at distance 
z [for example, Eq. (38) of Ref. [37], replacing rii with 
2Ndz, and then integrating over all planes z containing 
impurities. This procedure gives 



1 



Airhe 2 



2Ndz / dO 



(1 — cos ( 



(23) 

In this equation, kf = aKsfi/e 2 is the Fermi wave- 
length, <j>i(q;z) = (2ire 2 /K s q)exv{-qz}/{l + (qrs)- 1 } is 
the screened potential (in momentum space) created by 
a single impurity at position z, and q — 2fc/ sfn(#/2) is 
the change in momentum associated with scattering by 
an angle 0. 



Evaluating the integral of Eq. ( 23 ) at small a gives 



1 



ira In (1/a) 



2 N 



(24) 



Inserting this result for t into Eq. ( 22 ) and substituting 
/i = e 2 kf/aK s and kf 



Aim yields the result for con- 
ductivity announced at the beginning of the section, Eq 
pl. 



Equation ( 20 1 can be contrasted with the w idely-used 
result for the 2D model of charge impuritie d 36 * 37 * 39 ^, 
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for which the conductivity is linearly proportional to the 
electron density: a/(e 2 /h) ~ (l/a 2 )(n/ni). This dif- 
ference can be understood conceptually by noting that, 
for large angle scattering, only those impurities at a dis- 
tance smaller than the Fermi wavelength, Ay ~ n -1 / 2 , 
contribute significantly to scattering. One can therefore 
define, roughly speaking, an effective 2D concentration 
of scattering impurities as iVAy ~ N/n 1 / 2 . Inserting 
N/n 1 / 2 for rii gives a oc (l/a 2 )(n 3 / 2 /N), similar to Eq. 
( 20 ) . The remaining factor l/ln(l/a) in Eq. (20) is re- 
lated to low-angle scattering by distant impurities with 
z Xf. So far we are unaware of any transport data for 
TIs that shows a oc r?l 2 . Recent conductivity measure- 
ments on ultra-thin TIs (with thickness ~ 10 nm -C Ay) 
suggest 4 ^ a oc n, consistent with the 2D model of impu- 
rities. 

Our 3D model also produces a distinct result for the 
minimum conductivity <r m i n that appears in the limit of 
small average electron concentration. At small enough 
chemical potential that fi <C e 2 N 1 / 3 /a 2 ^ 3 K s , the sur- 
face breaks into electron and hole puddles, and one can 
think that the effective carrier concentration saturates at 
~ n p [see Eq. (10)]. An estimate of cr m i n can therefore 
be obtained by setting n ~ n p in Eq. (|20j), which gives^ 
the result of Eq. (21). 2D models oiaisorder impuri- 



ties also produce a minimum conductivity that is inde- 
pendent of the impurity concentration, but which has a 
different de pend ence on a. Specifically, at small a such 
models giv e 1 36 1 42 1 a min ~ (e 2 //i) ln(l/a). Our model sug- 
gests a minimum conductivity that is larger by a factor 
~ [aln^l/a)]- 1 . 



VIII. TI SURFACE WITH A GAP 



In Sees. \V\ and |VI[ we discussed the disorder potential 
created by Coulomb impurities at a gapless TI surface, 
whose massless spectrum is protected by time-reversal 
symmetry. On the other hand, one can open a gap at 
the TI surface if one introduces some source of time- 





FIG. 7. (Color online) Schematic illustration of a gap opening 
at the TI surface between the conduction band (upper, red) 
and valence band (lower, blue). With the addition of some 
source of time-reversal symmetry breaking, the ungapped dis- 
persion relation (left) acquires an energy gap U (right). 



is well-described by the results of the previous section. 
In particular, the disorder potential width r = Tq = 
(2 1 / 6 ^/a 2 / 3 )(e 2 N 1 ^/K s ) [see Eq. 0] and the average 
density of states (i/) = (v) = (a^ 3 /2^ 3 Tr)(K s N^ 3 /e 2 ) 



[see Eq. jl6[)]. If the gap U is small enough that U <C r , 
then the disorder potential at the surface is essentially 
unaffected by the gap, since local fluctuations in the 
Fermi level are much larger than U . For example, if the 
gapless surface spectrum is replaced with a "massive" 
dispersion relation, 



E = ± v / (hvk) 2 + {U/2f 



(25) 



as plotted in Fig. [7J then one can estimate the first-order 
effect of the gap by carrying out the same self-consistent 
procedure outlined in SeclVl In particular, the gapped 
dispersion relation of Eq. ( |25| has a corresponding den- 
sity of states 



(26) 



where @(x) is the Heaviside step function. At U /T <C 1, 
one can assume a Gaussian distribution of the Coulomb 



rive at a value for T [see Eq. (15)]. Expanding the result 
of this procedure for small U/Tq gives for the disorder 
potential width a slightly enhanced value 



reversal symmetry breaking, such as an external mag- potential <j> with some unknown variance T 2 , integrate 
netic fiel d 4 ' 43 ' 44 ', p roxim ity to a magnetic material or this distribution over <f> multiplied by u(-etf>; U) to pro- 
magnetic impurities 35 ™ the proximity effect from an ad- duce the thermodynamic density of states (v) , and then 
jacent superconducto ^, or electron tunneling between use the self-consistency relation T 2 = e 2 N/n s {v) to ar- 
two nearby TI surfaces^ 4 ^ (see also the review of Ref. 
[S]). The resulting gapped spectrum is illustrated schemat- 
ically in Fig. [7] 

In this final section we briefly discuss how the presence 
of a gap with magnitude U affects the disorder potential 
at the surface and the mid-gap density of states. We 
focus our discussion around the case where the chemi- 
cal potential fi = 0, which roughly corresponds to the 
largest disorder potential and the minimum in the ther- 
modynamic density of states. Again, we limit our con- 
sideration to the case where the TI is sufficiently thick 
that one can describe impurities as three-dimensionally 
distributed. 

In the absence of a gap, [7 = 0, the disorder potential 



r(U) ~ r i 



u 2 
24r 2 



(27) 



Similarly, the thermodynamic density of states is slightly 
depleted: 



U 2 



(28) 



On the other hand, if U is much larger than r , then 
the surface screens poorly and the disorder potential 
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grows. In this case screening of the disorder potential 
by the surface happens only nonlinearly, through the for- 
mation of electron and hole puddles at locations where 
the magnitude of the Coulomb potential energy reaches 
the gap energy U/2. This is similar to the bulk nonlin- 
ear screening discussed in Sec. |TlJ and naturally produces 
T(U) ~ U. The typical correlation length of the disorder 
potential at the surface (the nonlinear screening length) 
is given by 



Notice that if the gap U is reduced to the point where 
U ~ To, then (y) — > (v)o, as can be seen by comparing 
Eq. pi} with Eqs. || and 

Of course, these estimates assume that the surface gap 
U is smaller than the bulk band gap E g , and consequently 
that Rjj -C R, so that impurities near the surface are not 
screened by bending of the bulk bands. If the surface gap 
U is larger than E g , then the disorder potential variance 
is truncated at T(U) ~ Eg due to bulk screening. 

We note that Eqs. (301 and (31 1 were first derived 



Ri 



TVe 4 



(29) 



as in Eq. with E g -> U. 

One can estimate the corresponding concentration of 
electrons/holes in surface puddles, n p , by noting that 
a square area of size R\j at the surface should contain 
enough electrons/holes to neutralize the net charge of V — 2[i and 



Coulomb impurities in the adjacent cubic volume Rfj of 
the TI bulk. This gives n v R 2 v 
words 



in Ref. 50 in the context of semiconductor heterostruc- 
tures in a transverse magnetic field, where a 2D elec- 
tron gas experiences disorder from adjacent 3D impu- 
rities and the gap U in the kinetic energy spectrum 
is provided by the Landau level spacing fvjj c . These 
authors also showed how the disorder potential is re- 
duced and the density of states increased as the chem- 
ical potential fi is increased from zercP^. Specifically, 
e 2 N/n s {U ~ 2/i) 2 , provided that 



W NRfj, or in other 



e 2 N 

K„U' 



(30) 



The corresponding thermodynamic density of states 
can be estimated by noting that when the chemical po- 
tential [i is raised by an amount ~ U/2, surface hole 
puddles should dry up and be replaced by a correspond- 
ingly increased number of electron puddles. This suggests 
(y) = dfi/dn ~ n p /U, which gives 



e 2 N 
^U 2 ' 



(31) 



U - 2fi > e 2 N 1 / z /n s . Of course, for TIs the effect of a 
transverse magnetic field goe s bey ond simply opening a 
single gap at the Dirac poin#21Ml. We do not consider 
here the full problem of screening of Coulomb impurities 
in the presence of a magnetic field, but in principle this 
problem can be dealt with along the lines of Ref. l50l 
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